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ABSTRACT 

We present the validation and characterization of Kepler-61b: a 2.15 i?® planet orbiting near the 
inner edge of the habitable zone of a low-mass star. Our characterization of the host star Kepler-61 is 
based upon a comparison with the set of spectroscopically similar stars with directly-measured radii 
and temperatures. We apply a stellar prior drawn from the weighted mean of these properties, in 
tandem with the Kepler photometry, to infer a planetary radius for Kepler-61b of 2.15±0.13 R^ and 
an equilibrium temperature of 273±13 K (given its period of 59. 87756±0. 00020 days and assuming a 
planetary albedo of 0.3). The technique of leveraging the physical properties of nearby "proxy" stars 
allows for an independent check on stellar characterization via the traditional measurements with 
stellar spectra and evolutionary models. In this case, such a check had implications for the putative 
habitability of Kepler-61b; the planet is 10% warmer and larger than inferred from i^-band spectral 
characterization. From the Kepler photometry, we estimate a stellar rotation period of 36 days, which 
implies a stellar age of > 1 Gyr. We summarize the evidence for the planetary nature of the Kepler- 
61 transit signal, which we conclude is 30,000 times more likely to be due to a planet than a blend 
scenario. Finally, we discuss possible compositions for Kepler-61b with a comparison to theoretical 
models as well as to known exoplanets with similar radii and dynamically measured masses. 
Keywords: eclipses — stars: planetary systems — stars: individual (Kepler-61, KOI 1361, KIC 
6960913) 



1. INTRODUCTION 

With the discov eries of exoplanets Kepler-22b 
(iBorucki et all 120111) . Kepler- 20e fc f (iFressin et al.l 
[201J), Kep ler-42a, b, fee dMmriiead et all l2012aD. 
KeDler-68c (fGilliland et all I2013n. and KeDler-62e & f 
(jBorucki et al.l[2013l ). astronomers are encroaching upon 
the regime of transiting terrestrial exoplanets in their 
stellar habitable zones. Kepler-22b is the first super- 
Earth-sized exoplanet with a measured radius to reside 
in the habitable zone of a sun-like star, though its radius 
of 2.4 J?® does not necessitate a terrestrial composition. 
The Kepler-20, Kepler-42, and Kepler-68 exoplanetary 
systems each comprise multiple planets, some of which 
are Earth-sized or smaller (as small as Mars in the case 
of Kepler-42c). However, these planets orbit too close 
to their host stars to lie within the habitable zone. The 
star Kepler-62 hosts five planets, two of which are both 
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very likely solid and reside in their star's habitable zone 
(jBorucki et al.l[20T3l l. Th e most recent release of Kepler 
exoplanetary candidates (jBatalha et al.l 120131 ) contains 
10 members <2 R^ and with equilibrium temperatures 
between 185 and 303 K. This temperature range is a 
generous definition of the habitable zone proposed by 
iKastind (|2012( ). Half of these candidate exoplanets or- 
bit stars cooler than 4100 K, as reported by the Kepler 
Input Catalog. 

However, inferring the properties of low-mass stars 
from spectra (upon which a measurement of planetary 
radius and equilibrium temperature hinges so critically) 
presents difficulties on multiple fronts. The direct com- 
parison of theoretical spectra to observations, which is 
robust for deducing the properties of solar-type stars, 
is challenging for low-mass stars. Such spectra rely on 
detailed, computationally intensi ve modeling of convec- 
tion i n low-mass stell ar interiors (|Mullan fc MacDonaldl 
120011 : lBrowning|l2008l ) and complete lists of the complex 
array of molecules a nd grai ns that re side in their atmo- 
spheres (jTsuji et al . 1996; Allard et a l. 2000). For this 
reason, we often appeal to empirical, rather than theo- 
retical, method s for the phy sical characterization of low- 
mass stars fsee lTorreil2011l for a complete review). This 
challenge is compounded by the possibility that stellar 
properties may also depend on other parameters, such as 
activity and metallicity, in a significant way. The empir- 
ical technique for deriving M dwarf temperatures and 
m etallicities from K-hanA spectra that was innovated 
bv lRoias-Avala et al.l (J2012D offers an important inroad. 
However, this technique depends upon the H2O-K2 spec- 
tral index, which is effective probe of stellar temperature 
for mid-M dwarfs, but saturates for stars with tempera- 



tures higher than 3900 K (jMuirhead et al.ll2012bD . There 
exists a desert in stellar temperature, near 4000 K, where 
no reliable method method for temperature derivation 
from a spectrum exists: the H2O-K2 index has saturated, 
and the star is yet too cool for comparison between high- 
resolution optical spectra and synthetic models. This 
is an especially salient problem for the characterization 
of the current and future sample of exoplanets orbiting 
low-mass stars, given the astonishing occurrence rates 
of 1.0 planet/star (found from the Kepler sample to be 
O.90I003 planets/star per lDressing fc Charbonneaul2013l 
and l.OiO.l planets/star per lSwift et"aDl2013[ ). 

There exists an alternative means of measuring the 
properties of the nearest and brightest low-mass stars, 
with measured distances from parallax: interferometric 
measurements of their radii. The radius measurement, 
in tandem with the bolometric flux, also enables a di- 
rect measurement of the stellar temperature with min- 
imal modeling uncertainties. The number of low-mass 
stars with directly measured properties from interferom- 
etry is growing, and this pool of stars can be plumbed 
for proxies to sta rs too faint for such dire ct characteriza- 
tion themselves. iMuirhead et al.l (|2012al ) undertook the 
first such analysis with an application of the properties 
of Barnard's star toward a characterization of the M4V 
star and three transiting planets comprising the Kcplcr- 
42 system. A similar method would be especially useful 
for a star astride the 4000 K boundary, where spectro- 
scopic estimates of the stellar properties may be unre- 
liable, and for an exoplanetary host star. Kepler-61 is 
such a star: its temperature lies near to 4000 K, and de- 
pending upon the source of its stellar characterization, 
the equilibrium temperature of its planet lies either out- 
side or astride the stellar habitable zo ne. The radius 
of the plai iet varies from 2.0 to 2.3 i?® (jMuirhead et al.l 
l2012al and lBurke et al.l[20T3l respectively) depending on 
the assumed size of the star, which range brackets both 
a plausible rocky composition (more amenable to habit- 
ability) or a "mini-Neptune" composition. Moreover, the 
now exists a sample of four stars with (a) similar spec- 
tral type and (b) direct radius and temperature measure- 
ments, which can be applied to break the degeneracy of 
the planet's putative habitability. 

The apparent magnitude of Kepler-61, with Kepler 
magnitude Kp=15.0, renders the star too dim to en- 
able a radial velocity measurement of the planet's mass. 
However, even without a mass measurement of Kepler- 
61b, we are able to validate its authentic planetary 
nature with a statistical argument about the likeli- 
hood of the planet scenario in comparison to false- 
positive scenarios. We undertake such an analysis 
with BLENDER, which mach inery has already been ap- 
plied to validate Kepler-9d (jTorres et al.l 1201 ih. Keoler- 
llg jLissauer et al.l 
mm. KeDler -19b (j 



20111). Keple r -lOc J P Yessin et al.l 
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(iBorucki et"aIll2012D . and Kepler-20 e & f (jFressin et al.l 
I2OI2D . In this case, a single observation of the transit 
depth at 4.5 /im with Warm Spitzer in tandem with the 
BLENDER result, plays a prominent role in ruling out hi- 
erarchical triple false-positive scenarios. 

In Section 2, we describe the Kepler observations of 
Kepler-61. In Section 3 we describe our characterization 
of the transit light curve and the physical parameters of 



the star, including our method of applying the proper- 
ties of nearby similar stars to characterize Kepler-61. In 
Section 4, we describe the validation of Kepler-61b as an 
authentic planet with BLENDER. We include a description 
of follow-up imaging observations of the star to charac- 
terize any additional stars within the Kepler photometric 
aperture of Kepler-61, as measurements of the photocen- 
troid gather from the Kepler images and a measurement 
of the transit depth with Warm Spitzer. In Section 5, we 
comment on the stellar rotation of Kepler-61 and tran- 
sit times of Kepler-61b. We also discuss plausible com- 
positions for the planet, given the growing number of 
transiting super-Earth-sized planets with dynamically- 
measured masses. Finally, in Section 6, we describe fu- 
ture prospects. 

2. KEPLER OBSERVATIONS 

The Kepler spacecraft, launched on 7 March 2009, is 
photometrically monitoring 170 ,000 stars for 8 years for 
evidence of transiting planets. lArgabright et al.l ()2008l ) 
provides an overview of the Kepler instrume nt, and 
iCaldwenetaH (2010) and J enkins et al.l (|2010bl ) provide 
a summary of its performance since launch. The Kepler 
observations of Kepler-61 {Kepler Input Catalog num- 
ber 6960913) that we present in this work were gathered 
from 13 May 2009 to 3 October 2012, spanning Kepler 
"Quarters" 1-14. All data for this star were gathered in 
long-cadence mode (characterized by an exposure time 
of 29.5 minutes) for Quarters 1-11, and in short-cadence 
mode (characterized by an exposure time of 58.5 s) for 
Quarters 12, 13, and 14. The data contain gaps of ap- 
proximately 3 days between quarters for scheduled down- 
links. Kepler-61b was first identified as exoplanetary 
candidate Kepler O bject of Interest (KOI) 1361.01 by 
IBorucki et al.l ()2011|) . We employed the light curves gen- 
erated by the Kepler aperture phot o metry (PDC-Map) 
pipeline, described in lTwicken et al.l ()2010[ ). to which we 
add an additional step. We remove the effects of baseline 
drift by individually normalizing each transit as follows. 
We fit a linear function of time to the fiux immediately 
before and after transit (specifically, from 9 hours to 20 
minutes before first contact, equal to 2.5 transit dura- 
tions, and an equivalent time after fourth contact). 

3. ANALYSIS 

3.1. Derivation of Stellar Parameters 

While the physical characterization of isolated 
low-mass stars is a notoriously difficult problem 
(jSegransan et al.ir2003l : iTorresI 120111) , several recent tech- 
niques have offered promising inroads by tying spec- 
tra of M dwarfs to directly measured quantities. Such 
characterization of low-mass stars is crucial, as the 
Kepler mission has demonstrated that low-mass stars 
are hosts to sm all planets at a rate of 1.0 plane t/star 
(jDressing fc Charbonneau 2013: Swift et al.ll2013D . 

3.1.1. Characterization in the Literature 

Keple r-61 is classified as an MO star bv IMuirhead et aD 
(|2012b( ) , who employed iiT-band spectra of the star to in- 
fer stellar properties (in this work, they refer to Kepler- 
61b by its KOI notation, KOI 1361.01). They mea- 
sure ratios of equivalent widths of Na and Ca to de- 
termine the stellar metallicity, and deformation between 



continuum regions within j^-band (t he H2O-K2 index, 
fir st developed by iCovev et al.l 120071 and re-calibrated 
bv lRoias-Avala et al.ll2012D . which they interpolate onto 
the theoretical metallicity and H2 O-K2 surface from the 
BT-Settl late-type model spectra (jAllard et al.ll2012D to 
determine the stellar effecti ve temperature. T h e met al- 
licity relation published by iRoias-Avala et al.l ()2012|) is 
calibrated using b inaries comprising an M dwarf and an 
FGK star, which iMuirhead et al.l ()2012b[ ) then applied 
to a set of isolated low-mass stars in the Kepler sam- 
ple. For Kepler-61, they find Toff of 39291^^5 K and 
metallicity [Fe/H] of -0.02±0.11 (though they caution 
that applying the iiT-band metallicity metric to stars 
with temperatures higher than 3900 K rel ies upon an 
extrapolation of the iRojas-Ayala et al1l2012l metric) . By 
comparing this temperature and me tallicity to the Dart - 
mouth stellar evolutionary models (Po tter et al.l [20081 ) . 
they determine a stellar radius i?*= 0.55±0.07 Rq and a 
mass M^= 0.57±0.07 Rq. However, as described above, 
this location in temperature space is close to the loca- 
tion where the H2O-K2 index saturates (for Tg-g- > 3900 
K. IMuirhead et al]|2012b[ ). and the deformation between 
continuum region s is too small to effectively probe stel- 
lar temperature. iDressing fc Charbonneaul ()2013[ ) cor- 
roborates the result that temperatures derived from 
the H2O-K2 index are significantly lower than those 
derived from the comparison of the broadband colors 
to models for stars nea r the 3900 K marker. The 
IDressing fc Charbonneaul (|2013D method relies upon a 
comparison of the measured magnitudes of the star from 
the Kepter Input Catalog (2MASS JHK and Sloan filters 
g, r, i, and z) against the colors predicted from the Dart - 
mouth stellar evolutionary models (Pott er et al.l [20081 ). 
They assign a prior on stellar metallicity based upon the 
metallicity distribution o f the M dwarfs observed in the 
iCasagrande et ahl ()2008[ ) sample and a prior on height 
above the galactic midplane similar to that applied by 
iBrown e"Fa l. (2001) for the Kepler Input Catalog. For 
the colors of Kepler-61, they estimate an effective tem- 



perature of 4060^1^^^ K, a radius of 0.57;^;^'j' Rq, and 
a mass of 0.57 



0.11 
~^.q'qq Mq (C. Pressing, private communi- 



cation) . 

3.1.2. Optical and Near-Infrared Spectroscopy 

Using the FAST spectrograph on the 1.5 m telescope at 
Mount Hopkins, AZ, we gathered a spectrum of Kepler- 
61 in the range 5560-7570 A with 0.75 A resolution (we 
employed an integration time of 20 minutes to achieve 
a signal-to-noise ratio of 30). In Figure [31 we show 
the Kepler-61 spectrum in comparison to spectra of two 
nearby K7V stars, GJ 380 and GJ 820B, gathered with 
the same instrument (we observed the former on 24 April 
2012 with integration time 10 s for a signal-to-noise ra- 
tio of 210 pixel"^ and gathered the latter spectrum from 
the FAST Spectrograph Archive. It was observed on 4 
November 2010 at resolution of 1.5 A and an integration 
time of 2 s for a signal-to- noise ratio of 200 pixel" '^). 
We have denoted the wavelength regions employed by 
spectral typing sof tware, "The Hammer", developed by 
ICovev et al.l ()2007[ ). which we use to deterinine a spec- 
tral type of K7V for the KOI. iLeoine et al.l ([20TI) find 
that classification with the Hammer agrees with stellar 
classification from spectral indices within 1.0 subtypes. 



We consider the set of similarly typed stars with di- 
rectly measured radii. While spectral types for stars 
in the M0-K7 range compiled in the literature often 
vary by 1-2 subtypes, we def er to the spectral types 
listed in Boyaiian et al. (2012"). In this case, because 
the Pressing fc Charbonneau (2013) predicted temper- 
ature for Kepler-61 encompasses an effective tempera- 
ture as high as 4160 K (within la), we elect to ex- 
clude K5 stars from the sample of spectroscopically sim- 
ilar stars with resolved radii: the K5 stars GJ 820A 
and GJ 720B have effective temperatures of 4361ifcl7 
(IKervella fc Fouauel[200l [^n Belle fc von Brauhl 120091: 
Bovaiian et al.l I2012D and 4393±149 (jBovaiian et al.l 



2012f ) respectively, which lowest temperature estimates 



within Icr are still hotter than the la range predicted for 
temperature of Kepler-61 from broadband photometry. 
The sample of K7 and MO stars with radius measure- 
ments currently comprises four stars: GJ 380 (radius and 



temperatures measured gath ered by iLane et al. 



Ivan Belle fc von BraunI [20091 an dlBovaiian et alT 



GJ 338A (measurements fr om iBovaiian et al.l 



I2OOI 



201^ 
2012), 



GJ 338B (measurements from Boyaiian et al. '2012"), and 
GJ 820B (measurements from K crvclla & Fouquc 200^ 
and Ivan Belle fc von BraunI I2009D . We list the proper- 
ties of these stars in Table 1. 

We have compiled the set of i^-band spectra for these 
stars in addition to the KOI, which we depict in Fig- 
ure [TJ in order to estimate their [Fe /H] metallic i ties in 
a uniform fashion from the metric of IMann et al.l ()2013f ) 
(which metric is valid for stars as early as spectral type 
K5). Our spectrum for Ke p ler-61 w as gathered and pub- 
lished by IMuirhead et al.l ()2012aD with the Trip leSpcc 
instrument at Palomar Observatory (jHerter et al.| [2008'). 
They employed an exposure time of 6 minutes on 6 June 
2011, with resolution of iA {R of 2700), to obtain a 
signal-to- noise ratio of 60 pixel" ^. We gathered spec- 
tra for GJ 338A and GJ 338B on 27 Jan 2013 with the 
SpeX instr ument at NASA's Infrared Telescope Facil- 
ity (IRTF, iRavneret aL[[200l . and for GJ 380 on 17 
Pecember 2012. We employed the ShortXP observing 
mode (resolution of 5A, R of 2000) and exposure times 
of Is, Is, and 5s, respectively, for signal-to-noise ratios of 
150, 150, and 700 pixel"^. We gathered ou r spectrum of 
GJ 8 20B from the IRTF Spectral Library (jRavner et al.l 
I2009D . It was observed on 2001 October 20 with the 
same resolution, and has a signal-to-noise ratio of 800 
pixel" ^. We compute uncertainties on these metallicities 
from adding the sc atter in the J^-band metric quoted by 
IMann et al.l ()2013D of 0.11 dex in quadrature to the in- 
trinsic uncertainty in the value of [Fe/H]x from the error 
in the spectrum at the wavelengths that are operative for 
the metric. We list these derived metallicities in Table 
1, but note that no standard star possesses a metallic- 
ity consistent with la of the most probable [Fe/H] value 
for Kepler-61 of 0.03. Unlike in the case of Kepler-42 
(jMuirhead et al.ll2012b[ ). none of these nearby stars pos- 
sesses features consistent enough with the KOI to render 
one of them an single ideal "proxy" star for the Kepler 
target. 

3.1.3. Activity and Age Indicators 

We consider the activity levels of the stars in this sam- 
ple, in comparison to Kepler-61. We have measured the 



equivalent widths in Ha for both Kepler-61 from th e 
FAST spectra to those compiled by iGizis et"all ()2002[ ). 
They exhibit equivalent widths in Ha (all in absorption) 
ranging from -0.82 (for Kepler-61) to -0.50 (for GJ 338A). 
In addition, we have measured the rotation period for 
Kepler-61 from the Kepler photometry, which we com- 
pare to the rot ation per i ods li sted for three of the four 
nearby stars in iBarnesI (|2007f ) . The "gyrochronology" 
technique of mapping the rota tional period of a st ar to 
its age is described in detail in 'Barnesj (J20031 I2007D , and 
is used specifically in Barnes (2007) to estimate ages of 
1.96 Myr, 1.36 Myr, and 2.96 Gyr for GJ 380, GJ 338A, 
and GJ 820, respectively. In Figure [U we depict the first 
four quarters of Kepler observations of Kepler-61, or ap- 
proximately one year of continuous observation. These 
observations have been processed using the Presearch 
Data Conditioning (PDC) module of the Kepler data 
analysis pipeline, with the Bayesian Maxim um A Poste- 
riori (MAP) approach applied (described in lSmith et al.l 
120121) . The use of highly correlated and quiet stars to cre- 
ate a set of co-trending basis vectors enables the removal 
of non-astrophysical artifacts from the Kepler time series, 
and the preservation of astrophysically interesting sig- 
nals such as stellar r otation. We apply the d iscrete cor- 
relation fu nction of lEdelson fc KrolikI (|1988[ ). similarly 
applied by iFabrvckv et al.l ( 20121) ori the time series of 
Kepler-30 and lQueloz et all (|2009[ ) on CoRoT-7, on this 
portion of the Kepler-61 light curve. We test lags from 1 
to 100 days, and identify a 36±4 day periodicity, which 
we attribute to the stellar rotation period. In Figure [21 
we show both the Kepler photometry and autocorrela- 
tion funct ion that we employ ed to characterize the stellar 
rotation. Ilrwin et al.l ()2011f) recently published a com- 
pilation of the known rotation periods of low-mass stars 
from the literature. The rotation periods are drawn from 
open clusters of stars with derived ages from 1-650 Myr 
(which ages are measured by main sequence fitting of 
these clusters, as compared to stellar evolutionary mod- 
els), and then also from field stars with ages >1 Gyr. 
Among the clusters with ages less than 1 Gyr, nearly all 
of the stars with masses >0.5 Mq have rotation periods 
shorter than 30 days, and only after 1 Gyr do low-mass 
stars with masses greater than 0.5 Mq appear to spin 
down enough to produce 36 day rotation periods. Stars 
with masses between 0.5 and 0.7 Mq and ages between 8 
and 10 Gyr are more likely to have rotation periods in the 
tens of days (Kiraga & Stcpicn 2007 and Baliuna s et al.l 
119961 observed values <30 days for field stars in the 1- 
2 Gyr range). We therefore take the observed rotation 
period of Kepler-61 to be conservatively indicative of an 
age >1 Gyr. 

We also apply the age metric of IBarnesI (I2010f) to esti- 
mate the age of Kepler-61, which lSwift et al.l ()2013f ) used 
to obtain an age approximation for Kepler-32: 
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where the dimensionless constants kc = 0.646 da ys Myr^ 
and k f = 452 Myr day^ are approximated in IBarnesI 
(|201CI|). We estimate t he convective turnover time r 
from I Wright et ahl ()2011l ) , who derived empirical convec- 
tive turnover times for a sample of 824 solar and late- 
type stars, the typical convective turnover time of a star 



within the mass range 0.47-0.62 Mq is 29 days. It's 
not possible to determine the initial spin period Pq as- 
sociated with Kepler-61 at its birth, so we defer to the 
median value of the initial spin period, Pq = 2.81 days 
required to produce the observed rotation rates for 0.6 
Mq stars i n the Praes e pe clu ster (jAgiieros et al.l [20Tll ) 
similarly to iSwift et al.l ()2013[ ) . This formulation returns 
an age for Kepler-61 of 10 Gyr. We therefore find that 
ages between 1-10 Gyr are consistent with different met- 
rics of age constraints for Kepler-61, and simply adopt a 
lower bound on its age of 1 Gyr. 

Because none of these nearby stars comprises an ideal 
"proxy" to the planet-host star (i.e., possessing both sta- 
tistically indistinguishable metallicities and temperature 
indices in tandem with similar activity indicators), we 
adopt the conservative tack of employing a radius and 
temperature for Kepler-61 which are the weighted mean 
of the radii and temperatures of the set of standard stars. 
For our uncertainty on these values, we encompass the 
highest and lowest mean value among the sample. We 
therefore adopt for Kepler-61 a radius of 0.62lg'Q5 Rq 
and a temperature of 4017 l^^gQ K. 

3.2. Derivation of Planetary Parameters 

We estimated the uncertainty in the planetary tran- 
sit parameters using the Markov Chain Monte Carlo 
(MCMC) method as follows. We emp loy model light 
curves generated with the routines in IMandel fc Agoll 
(2002), which depend upon the period P, the epoch Tc, 
the planet-to-star radius ratio Rp/Ri,, the ratio of the 
semi-major axis to the stellar radius a/R^,, the impact 
parameter &, the eccentricity e, and the longitude of pe- 
riastron, uj. We fixed two quadratic limb-darkening co- 
efficients (LDCs), Ml and U2, to theoretical values based 
on the adopted effective temperature and radius. We 
employed the theoretical limb darkeni ng coefficients gen- 
erated for the Kepler bandpass by iClaret fc Bloemeiil 
(|2011f) from the PHOENIX models corresponding to a 
star with Tes of 4000 K and solar metallicity, which is 
the closest match to the effective stellar temperature Teg 
of 4017 and the metalhcity of 0.03±0.14 for Kepler-61 
(we additionally specified a log(g) of 4.5 and an interme- 
diate turbulent velocity value of 2 km s~^): these coeffi- 
cients are ui=0.50 U2=0.20. We accounted for the 29.5 
minute integration time of the Kepler photometry, which 
is three times longer than the ingress and egress duration 
of the planet candidate, by evaluating the light curve 
model at intervals of 1 minute, and then summing the 
model over the long cadence integration time. We model 
the three quarters of short cadence observations inde- 
pendently (gathered from Quarters 12-14, during which 
time the planet presented five transits). For these short 
cadence observations, we evaluate the light curve model 
at each time measurement, in 58.5 s intervals, and fix the 
period and transit time to the best- value recovered over 
the 11-quarter baseline of long-cadence observations. 

To generate the MCMC chain, we randomly choose 
one parameter, perturb it, and evaluate the x^ of the 
solution. If the x^ is lower, we accept the new parame- 
ter value. If the x^ is higher, we evaluate the probabil- 
ity of accepting the jump as p = e~ -^ '^. If the jump 
is rejected, the procedure is repeated at that point in 
the chain until an acceptable jump occurs. We adjust 



Table 1 

Comparison of Observables between Kepler-61 and Similar Stars 



Parameter 



GJ 380 



GJ 338A 



GJ 338B 



GJ 820B 



Kepler-61 



Spectral I'ype'' 

Metallicity [Fe/H] 

H2O-K2 

Ejj ^ (in absorption) 

Rotation period [days]'^ 
Estimate Age [Gyr]"^ 


K7V 

0.18±0.11 

1.044±0.002 

-0.61 

-5.16 

11.67 

0.2 


MOV 

-0.15±0.12 

1.034±0.005 

-0.56 

-4.68 

0.1 


Krv 

-0.15±0.12 

1.025±0.005 

-0.50 

-4.65 

10.17 


K7V 

-0.23±0.11 

1.019±0.002 

-0.59 

-5.03 

37.9/48 

3 


K7V 

0.03±0.14 

1.02 ±0.010 

-0.82 

36 



0.62+^5 (adopted) 
4017^^50 (adopted) 



R* [Rq]" 



0.64±0.004 
4085±14 



0.58±0.01 
3907±35 



0.57±0.01 
3867±35 



0.60±0.007 
3932±25 



^ Values for standard stars from lBovaiian et al.l I I2012I) . 

Values for standard stars from lGizis et all i2002r . 
'^ From lBarnesI J2007I) . stated to nearest 100 Myr. 

the width of the distribution from which we randomly 
draw the jump sizes in each parameter until 20-40% of 
jumps are executed in each of the parameters. We cre- 
ated five chains, each of length 10^ points, where each of 
the chains is begun from a different set of starting pa- 
rameters (each parameter is assigned a starting position 
that is +3a or -3a from the best-fit values). We dis- 
card the first 20% of jumps from each chain to remove 
the transient dependence of the chain on the starting pa- 
rameters. We first conducted this analysis, as described, 
using only the Kepler light curve to inform our value 
of x^- However, the allowable stellar densities we in- 
fer from the light curve alone are much broader than 
the range of stellar densities consistent with our spec- 
troscopic information about the star. Because Kepler- 
61 is a late K star, it is slowly evolving- therefore, its 
range of theoretical densities is tightly constrained for 
ages < 14 Gyr. When we apply the MCMC algorithm 
toward fitting the long-cadence transit parameters inde- 
pendently and allow a/Ri, to float, we find that values 
of a/Ri, from 35-150 furnish comparable fits to the light 
curve (the ingress and egress time, at 9.3±3.2 minutes 
measured at long-cadence, results in a wide family of 
allowable light-curve fits) . We take advantage of two cir- 
cumstances that allow us to better constrain the transit 
parameters. First, there exist 5 transits gathered by Ke- 
pler at short-cadence, where the ingress and egress time 
are resolved by the 1-minute exposure time. Second, we 
make use of the fact that low-mass stars are slowly evolv- 
ing to set a physically-motivated prior on a/R^, as follows. 
We based our procedure for constraining the mass, ra- 
di us, and age o f the h ost star on the method described 
bv [Torres et al.l ()2008() . Using the metallicity we derive 
from K band, we created a set of stellar ev olution mod- 
els fr om the Dartmouth isochrone series (' Potter et al.l 
I2008f ). We employed the interpolation software that ac- 
companied that work, which accepts as inputs the age 
of the star, the iron abundance, and the abundance of 
a-elements relative to solar (for which we assume the 
solar value) , and outputs a grid of stellar isochrones cor- 
responding to a range of masses. We evaluated a set of 
isochrones over an age range of 1 to 14 Gyr (at intervals 
of 0.1 Gyr) and in [Fe/H] in increments of 0.01 from -0.5- 
0.5 dex (encompassing 3a above and below the measured 
[Fe/H] of 0.03 ± 0.14). We evaluate the physical radius 
corresponding to each stellar model via log(g) and the 



mass of the star {g = GM^^/Rl). 

Rearranging Kepler's versio n of Newto n's thir d law 
in the manner employed b y iSeager fc^ allcn-O rnelasI 
(|2003f) . [Soizetti et al.l((2007D and lTorres et al.i (|2008l ). we 
convert the period (derived from photometry), and the 
radius and mass of the host star (from isochrones) to a 
ratio of the semi-major axis to the radius of the host star, 
a/Ri,: 
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(2) 



where we will hereafter assume that Mp is negligible 
when compared to the mass of the host star. 

We calculate the corresponding value of a/R. for each 
stellar model. We then generate the MCMC chain as 
follows. We implement a prior on a/ R^, by varying the 
adopted quantities i?* and T„a in the chain, in addition 
to the light curve parameters. For each set of i?* and T^a, 
we locate the closest stellar model associated with these 
values and record its corresponding stellar density, a/ R^,. 
It is this value of a/Ri, that is used to generate the light 
curve model, along with the other light curve parameters 
P, Tc, Rp/Ri,, b, e, and uj (defined above), which are 
permitted to vary independently. We assign uniform flat 
priors in b (from 0-1), e (from 0-1), and uj (from 0-27r), 
and uniform improper priors on all other parameters. In 
this way, we are sampling only values of a/ R^, that are 
consistent with the spectroscopically-derived parameters, 
but values of a/i?+ that are not as well matched to the 
light curve are penalized by the x^ term corresponding 
to the photometry. We adopt Gaussian priors on R^, 
and T„Q, which we implement by adding extra terms in 

the x^ (where P corresponds to the vector of light curve 
parameters at each iteration): 
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In Figure 13.21 we show the correlations between the 
posterior distributions of subset of parameters in the 
model fit, as well as the histograms corresponding to 
each parameter. In Figure |B1 we show the phased Ke- 
pler transit light curve for Kepler-61b, with the best-fit 
transit light curve overplotted. We report the best-fit 
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parameters and uncertainties in Table 2. The range of 
acceptable solutions for each of the light c urve param- 
eters is determined as follows. Following i Torres et alj 
(|2008f l. we report the most likely value from the mode 
of the posterior distribution, marginalizing over all other 
parameters. The uncertainty is derived from the extent 
of the posterior distribution that encloses 68% of values 
closest to the mode. 

3.3. Physical Parameters 

This procedure described in Section 1.1 is also advan- 
tageous in that, in addition to recording the value of 
a/Ri, at each iteration of the MCMC chain, we may also 
record the other traits of the star at that value of i?* 
and Tgff, including its mass, luminosity, and age. At the 
conclusion of the MCMC analysis, therefore, we have ac- 
cumulated a chain not only for the light curve parame- 
ters, but for the physical parameters as well, as predicted 
from evolutionary models. The correlations between pa- 
rameters, both physical parameters and those associated 
with the light curve, are therefore preserved in the chain 
and incorporated into our estimate of the stellar param- 
eters (although we note that we have not accounted for 
possible correlated error between our adopted values of 
effective temperature and radius measurements of the 
star). We calculate the planetary radius from multiply- 
ing the elements of the Rp/Ri, by the chain of i?^. We 
infer a value for the stellar mass from its posterior distri- 
bution of M*=0.635±0.037. The slowly-evolving nature 
of Kepler-61 results in a largely unconstrained estimate 
of stellar age. The stellar rotation period (36 days, de- 
scribed in Section 6.1) indicates a star older than 1 Gyr, 
which is consistent, though also itself only a weak con- 
straint. It's also possible to evaluate the posterior distri- 
bution of planetary equilibrium temperatures from the 
MCMC analysis. In the case of a circular orbit, we re- 
quire only the stellar radius and temperature, the plane- 
tary semimajor axis, and the planetary albedo. How- 
ever, in the case of an eccentric orbit, the planet re- 
ceives time-variable stellar insolation. In order to evalu- 
ate the equilibrium temperature of the planet in the case 
of non-zero eccentricity, we evaluate the time-averaged 
equilibrium temperature by performing an integral over 
the mean anomaly from to 27r, u sing the formalism 
detailed in lMurrav fc Correial (|201Cl[ ). where A is the ge- 
ometric albedo of the planet, d is its distance from the 
star, and M is the mean anomaly: 



U^-Af/^ r sI^TAM (4) 
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In Figure [51 we show both the planetary temperatures 
based on the apastron and periastron isolation from the 
star, as well as the time-averaged temperature of the 
planet, for each element of the MCMC chain. 

We imposed a fiat prior on the eccentricity of Keplcr-61 
from 0-1, so the posterior distribution on e includes val- 
ues as high as 0.9 (such a large eccentricity still matches 
the transit duration at finely-tuned values of w). Indeed, 
the circularization timescale for Kepler-61 has not yet 
elapsed if we assume it did not possess a large initial ec- 
centricity. We consider the expressio n for circularization 
time (for modest initial e) given by iGoldreich fc Soter, 



(|1966f) . where a is the semimajor axis of the planet, Rp 
is the planetary radius, Mp is the planetary mass, M^, 
the stellar mass, Q is the tidal quality factor for the 
planet (which is highly uncertain, but we test two values: 
100 for the assumption of a terrestrial composition, per 
Goldreich fc Soter 1966, or 10^, which is the lower limit 
measured bv lBanfield fc MurravllT99^ for Neptune) and 
G is the gravitational constant: 
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We find that the circularization timescale would be 400 
Myr for a terre strial Q of 100 (assuming a 7 M^ planet 
orbiting a 0.64 M© star at 0.25 AU), and 190 Gyr for a 
Neptune-like Q of 10^. It's therefore plausible that the 
planet resides in an eccentric orbit. If Kepler-61 began 
with a large initial eccentricity, then terms of order (1- 
e^) become relevant and the circularizatio n timescale de- 
creases (as elucidated in Equations 7-9 of iSocrates et al.l 
2012). It's therefore also possible that the circulariza- 
tion timescale has indeed elapsed, dependent upon the 
initial eccentricity and tidal Q of the planet. However, 
the uncertainty of our knowledge about its initial eccen- 
tricity (coupled with the uncertainty about the correct 
value of tidal Q for Kepler-61 in particular) is such that 
we believed a flat prior on e to be appropriate. 

We find consistent orbital parameters for the 11 quar- 
ters of long-cadence observations (13 transits) and the 
three quarters of short cadence observations (5 transits), 
though our ability to resolve the shape of ingress and 
egress with short-cadence mode, coupled with the prior 
on a/ Ri, from our knowledge of the stellar parameters, is 
reflected in the error bars on a/Ri, and the ingress/egress 
time T, even though we have only half the number of 
transits in the latter mode. While we find r=9.3±3.2 
minutes from the 13 transits in long-cadence mode, we 
find r=9.56±0.47 minutes from the five transits in short- 
cadence mode. 

4. PLANETARY VALIDATION OF KEPLER-61 



IMorton fc Johnsonl (j2011f ) provide a priori false pos- 
itive probabi lities for the Kepler planetary candidates 
published bv lBorucki et al.l (|2011[ ). within which sample 
Kepler-61 is included. They cite the vet ting of candi- 
dates by the Kepler software (detailed by iBatalha et al.l 
l2010f ) as being already sufficient to produce a robust list 
of candidates, and combine stellar population synthesis 
and galactic structure models to demonstrate that nearly 
all of these 1235 candidates have a false positive prob- 
ability <10%. Kepler-61b, with a Kepler magnitude of 
14.995 and a galactic latitude of 9.6°, has an a priori 
false positive probability of 4.8%. 

4.1. Adaptive Optics Imaging 

We place limits on the presence of additional stars in 
the neighborhood of Kepler-61 with adaptive optics (AO) 
observations gathered at Keck with the NIRC2 instru- 
ment. On 22 June 2012 we observed Kepler-61 in both 
J and K band, with observed FWHM of the core in K 
of 0.06" and in J of 0.10". We detect a companion star, 
2.9 magnitudes fainter and located 0.94" to the north- 
west of Kepler-61, which is shown in Figure [9] The high 



Table 2 

Star and Planet Parameters for Kepler-61 



Parameter 



Value & la confidence interval 



Keplcr-61 [star] 



Right ascension" 


19h41ml3.09s 


Declination'' 


-|-42d28m31.0s 


^cff [KP 


mnttto 


R* [Solar radii] >= 


n fi9-l"002 


M* [Solar masses 


0.635±0.037 


[Fe/H]^ 


0.03±0.14 


Age [Gyr 


>1 


Kepler-61 [planet 


Period [Days 


59.87756±0.00020 


To [BJD-2450000 


4984.1880iE;;'^S^! 


Rp/R-M 


0.03301±0.00085 


a/i?* 


90.6±3.4 


inc [deg 


>89.80 


e 


<0.25 


e cos(a;) 


0.0±0.29 


Impact Parameter 


< 0.29 


Total Duration [min 


290.7±4.6 


Ingress Duration [min 


9.56±0.47 


Rp [Earth radii 


2.15±0.13 


Planetary T^q [K 


273±13 



'^ ICRS (J2000) c oordinates from th e 2MASS AU-Sky Catalog 
of Point Sources iCutri et al.|[2003l) . The proper motion de- 
rived bv lRoeser et ah! I I2010I) is -5.6 milliarcscc yr^ in right 
ascension and 11.8 milliarcsec yr^ in declination (both with 
error bar of 3.8 milliarcsec yr~ ). 

Stellar temperature and radius inferred from weighted mean 
of directly measured K7 and MO sample, as described in text. 
^ Metallicity derived from ii"— band using metric from 
IMann et all I I2013I '). 



resolution i^T-band AO image cleanly resolves these two 
sources, while we employed the J-band image to bet- 
ter characterize the neighboring source. The additional 
source falls within the Kepler aperture, but is removed 
enough from the target star that we can employ the Ke- 
pler centroids to assert that the planet orbits the brighter 
star, as we describe in Section l473l The blended magni- 
tudes are 15.0 in the Kepler bandpass, J=13.077±0.022 
magnitudes and K— 12.272±0.019 magnitudes. We in- 
dependently measure the J and K magnitudes for the 
two stars, and apply the Kp-X and J-K relationships 
(which are derived separately for dwarfs and giants from 
the Kepler Input Catal og: this conversion i s described in 
detail in Appendix A of lHowell et al.|[2Q12D to determine 
the de-blended Kepler magnitude. For the primary tar- 
get, we find Kp = 15.22±0.09 mag, J = 13.149±0.022 
magnitudes, and K = 12.345±0.019 mag. For the sec- 
ondary, we find Kp= 18.20±0.10 mag, J=16.064±0.025 
mag, and /'ir=15.242±0.020 mag. The two stars exhibit 
indistinguishable J-K colors, with 0.804 ± 0.029 for 
Kepler-61 and 0.822±0.32 for the dimmer companion. 

We assess our sensitivity to additional source s using 
a simi lar procedure to that described by Batalh a et al.l 
(|201lD . We inject fake sources near the target star at ran- 
dom position angles, using steps in magnitude of 0.5 mag 
and varying the distance from the target star in incre- 
ments olF 1.0 FWHIVI of the point-spread function (PSF). 
We then attempt to id entify the inje cted sources with 
the DAOPhot routine fSt etsonI |1987|) and also by eye, 
and set our sensitivity limit, as a function of distance, at 
the magnitude where we are able to recover the injected 



sources. The limit in Am as a function of distance from 
the target star is shown in Figure [9l We then convert 
the Am sensitivity limit in K band to a limit in Kepler 
m agnitudes, again usi ng the Kp-K relationship detailed 
in lHowellet aLl(|20Tl ). 

4.2. Speckle Imaging 

We gathered speckle imaging of Kepler-61 on 11 June 
2011 UT using the two-color Differential Speckle Sur- 
vey Instrument at the Wisconsin Indiana Yale NOAO 
(WIYN ) 3.5 m telescope , located at Kitt Peak Obser- 
vatory (jHorch et al.ll2009l) . The speckle camera obtained 
7000 40 msec images in / band (8880/400 A). We reduced 
and processed these observations to produce a final re- 
constructed speckle image for the star. Details of the 
speckle camera observations for the Kepler follow-up ob- 
serving prog ram, including the reduction methods, are 
presented in lHowell et al.l ()2011[ ). 

The speckle observations allow detection of a compan- 
ion star within the approximately 2.76 x 2.76 arcsec box 
centered on the target. We can detect, or rule out, com- 
panions between 0.05 arcsec and 1.8 arcsec from Kepler- 
61 and, in this case, we did not detect a companion star. 
We report the limiting difference in magnitude for an 
additional star that would have been detectable with 3cr 
confidence in Table 3. The companion detected with 
adaptive optics imaging, which we describe in the pre- 
vious section, lies just beyond detectability with speckle 
imaging in / band, at a distance of 0.94" and Am in K 
band=2.9 (0.25 magnitudes from what would have been 
detected in the speckle image with 3(7 confidence at that 



distance from the star). 

4.3. Photocenter Tests 

We use two methods to search for false positives due to 
background ecUpsing binaries, based on examination of 
the pixels in the aperture of Kepler-61: direct measure- 
ment of the source location via difference images, and 
inference of the source location from photocenter motion 
associated with the transits. We employ two methods be- 
cause of their different vulnerabilities to systematic bias; 
when the methods agree, we have increased confidence 
in their result. 

Difference image analysis (jTorres et al.ll201ll) takes the 
difference between average in-transit pixel images and 
average out-of-transit i mages. A fit of the Kepler pixel 
response function (PRF; iBrvson et al.ll2010l ) to both the 
difference and out-of-transit images directly provides the 
location of the transit signal relative to the host star. 
We measure difference images separately in each quarter, 
and estimate the transit source location as the robust 
uncertainty-weighted average of the quarterly results. 

We measure photocenter motion by computing the 
flux- weighted centroid of the pixels in the optimal aper- 
ture, plus a one-pixel halo in every cadence, generating a 
centroid time series for row and column. We fit the mod- 
eled transit to the whitened centroid time series trans- 
formed into sky coordinates. We perform a single fit for 
all quarters, and then infer the source location by scaling 
the difference of the se two centroids by th e inverse of the 
flux as described in lJenkins et al.l ()2010al ). 

The source as determined by the difference image 
method is offset from the nominal location of Kepler- 
61, as given in the Kepler Input Catalog, by 0.09 ± 0.29 
arcsec = 0.68o-. The source as determined by the flux- 
weighted centroid method is offset from Kepler-61 by 0.32 
± 0.37 arcsec = 0.86cr. Both methods show that the 
observed centroid location is consistent with the transit 
occurring at the location of Kepler-61, and rule out the 
companion in the adaptive optics imaging as the source 
of the transit, which is 3cr removed from position at which 
the transit occurs. 

4.4. Spitzer Observations 

Warm Spitzer observations in the near-infrared can 
also prove useful toward validating Kepler candidates, 
as shown for Kepler-lO c (Fressin et al.;, 2 011i). Kepler- 14b 
(iBuchhaveet al.ll201lD. Kepler-18c fc d ijCochran et al.l 
120111). KeDi en9b (IBa llard et al.l [201l|), Kepler-22b 
(iBorucki et"all 1201 1[ ). Kepler-25b fc c (iSteffen et al.l 
20121). and K epler-20c. d, e, & f (IGautier et al.l 120121: 



Fressin et al.l 12012') Unless a putative blend scenario is 



comprised of stars of nearly identical color, the transit 
depth in a blend scenario will depend upon the wave- 
length at which it is observed. Conversely, an authentic 
transiting planet will produce an near-achromatic transit 
depth. 

We gathered observations using the Infrared Array 
Camera (IRAC) (Fazio et al,ii20Q4i) on Warm Spitzer at 
4.5 nm of the UT 17 September 2011 transit of Kepler- 
61b. The observations span 10 hours, centered on the 
4.75-hour-long transit. We gathered the observations us- 
ing the full-array mode of IRAC, with an integration time 
of 12 s/image. We employed the techniques described 



in lAgol et all ()2010l ) for the treatment of the images be- 
fore photometry. We first converted the Basic Calibrated 
Data products from the Spitzer IRAC pipeline (which ap- 
plies corrections for dark current, flat field variations, and 
detector non-linearity) from mega-Janskys per steradian 
to data number per second, using 0.1469 MJy-sr~^ per 
DN s~^, and then to electrons per second, using the gain 
of 3.71 e DN~^. We identified cosmic rays by perform- 
ing a pixel-by-pixel median filter, using a window of 10 
frames. We replace pixels that are > 4cr outliers within 
this window with the running median value. We also cor- 
rected for a striping artifact in the Warm Spitzer images, 
which occurred in the same set of columns, by taking the 
median of the pixel values in the affected columns (us- 
ing only rows without an overlying star) and normalizing 
this value to the median value of neighboring columns. 

We discuss several means of performing the Warm 
Spitzer IRAC photom e tric re duction for similar observa- 
tions in IBallard et al.l (|2011h . and make use of the con- 
clusions from that work. First, we estimate the position 
of the star on the array with a flux-weighted sum of the 
position within a circular aperture of 3 pixels. We then 
performed aperture photometry on the images using the 
centroid positions and variable aperture sizes between 
2.1 and 4.0 pixels, in increments of 0.1 pixels up to 2.7 
pixels, and then at 3.0 and 4.0 pixels. We decided to use 
the position estimates using a flux-weighted sum at an 
aperture of 2.6 pixels, which minimized the out-of-transit 
RMS. 

We remove the effect of the IRAC intrapixel sen- 
sitivity variations, or th e "pixel-phase" effect (see eg. 
iCharbonneau et all 120051 : iKnutson et al.l 120081 ) using a 
polynomial functional form for the intrapixel sensitivity 
(which depends upon the x and y position of the star on 
the array). We denote the transit light curve / (which 
depends upon time), and we hold all parameters constant 
except f or the transit depth. W e use the light curve soft- 
ware of iMandel &: Agoll (|2002l ) to generate the transit 
models. The model for the measured brightness f'{x,y) 
is given by: 



/' = f{t,Rp/R,y[bi+b2{x^x)+b3ix^xf+b^{y~y)+b5{y~yy 

(6) 

where we include all of the observations (both in- and 
out-of-transit) to fit the polynomial coefficients and the 
transit depth simultaneously. 

We fit for the polynomial coefficients 5i through 65 us- 
ing a Levenberg-Marquardt x^ minimization. However, 
the Spitzer light curve contains significant correlated 
noise even after the best intrapixel sensitivity model is 
removed. We incorporate the effect of remaining corre- 
lated noise with a resi dual permutation a nalysis of the 
errors as described by I Winn et"al] ()2008n , wherein we 
find the best-fit model /' to the light curve as given by 
Equation [51 subtract this model from the light curve, 
shift the residuals by one data point in time, add the 
same model back to the residuals, and refit the depth 
and pixel sensitivity coefficients. We wrap residuals from 
the end of the light curve to the beginning, and in this 
way we cycle through every permutation of the data. We 
determine the best value from the median of this distri- 
bution, and estimate the error from the closest 68% of 
values to the median. Using the residual permutation 



Table 3 

Magnitude Limits on Companions to Kcpler-61 from Speckle 
Imaging 

Radius of Annulus around Keplor-61 Limiting Delta Magnitude 
[Arcseconds] 3cr Confidence 



0.05-0.30 
0.30-0.50 
0.50-0.70 
0.70-0.90 
0.90-1.10 
1.10-1.30 
1.30-1.50 
1.50-1.70 
1.70-1.90 



2.69 
3.05 
3.09 
3.16 
3.15 
3.11 
3.18 
3.24 
3.20 



method on the hght curve treated with a polynomial, we 
find i?p/i?^,=:0.0315±0.0069, in excellent agreement with 
the Kepler measurement of Rp/Rstar=0.0M76±0. 00094. 
We note that t he use of the weight ed sensitivity func- 
tion proposed in IBallard et al.l (J2010D made a negligible 
difference to the photometric residuals in this case, so for 
reasons of computational time, we deferred to the poly- 
nomial reduction technique. In Figure 1101 we show the 
combined and binned Spitzerlight curve, with the best-fit 
transit model derived from the Spitzer observations and 
the best-fit Kepler transit model (with the quadratic limb 
darke ning coefficients for the S pitzer 4.5 fiia filter, drawn 
from IClaret fc BloemenI 1201 II as similarly described in 
Section 1) overplotted. 

4.5. BLENDER Validation 

Traditional confirmation of the planetary nature of a 
transit signal relied upon a dynamical mass measurement 
from radial velocity observations. In cases where the re- 
flex motion induced on the host star by the planet is too 
small to be detected, dynamical confirmation may yet be 
possible via transit timing variations (TTVs). However, 
in cases where a dynamical mass measurement via either 
technique is not viable, it's still possible to "validate" 
the planetary nature of the transit signal, via a statisti- 
cal argument about the relative likelihood of an authen- 
tic planet producing the transit signal, as compared to 
a false positive scenario. BLENDER is such a machinery, 
which combines evidence from the Kepler photometry 
(as compared to model light curves for planetary and 
false positive scenarios), spectroscopy, Spitzer photom- 
etry (where available), the stellar colors, and adaptive 
optics to deduce a false positive probability for a plan- 
etary candidate. BLENDER has already been applied to 
validate planets in a numbe r of Kepler exoplan etary sys- 
tems, including CoRoT-7 (iFressin et al.ll2012D. Kepler - 
lO lFressin e t al. (2 0111). Keple r-18 (iCoch ran et al.ll2011[) 
KeDler-19 (iBallard et al.ll2011h. K epler-20 (iGau tier et al.l 
[201^. KeDl er- 21 (IHowell et al.l [Ml . and Kepler-22 
(jBorucki et al.l I2Q12D . and its details are described 



therei n as w ell as in iTorres et al.l (|2Q04f ). IFressin et alj 
([2011'), and 'Fress in et al.l ( 20121) . We summarize the 
BLENDER procedure below, and present the statistical 
likelihood that the transit signal presented by Kepler- 
61b is attributable to a 2.5 i?® planet orbiting a 0.65 Rq 
star. 

The exquisite precision of the Kepler photometry is 
already sufficient to rule out some false-positive scenar- 



ios, which would produce a significantly different transit 
shape from the one observed by Kepler. Such a false 
positive could mimic the observed transit depth if an ad- 
ditional star fell within the same aperture of the Kepler 
target star. The light contributed by this undetected 
companion (which may be gravitationally bound to the 
target star or lie in the foreground or background) , would 
reduce the transit depth produced by an eclipsing binary 
system or a planetary system comprising a star and a 
larger planet, conspiring to produce a planetary tran- 
sit depth. BLENDER manufactures synthetic light curves 
corresponding to these false positive scenarios: for those 
with a physically bound companion, BLENDER assumes a 
common age for the putative companion star and the 
Kepler target star, while an unassociated background 
or foreground star is assigned an age of 3 Gyr. The 
mass of this secondary star and the tertiary body (either 
star or planet) is allowed to vary. These model blend 
light curves are compared to the Kepler photometry in 
a x^ sense. Blend scenarios that furnish a good fit to 
the Kepler light curve (within 3a of the best authen- 
tic planet model) are then tested for consistency against 
other constraints. These include (a) the color of the star 
as reported in the KIC (Brown et al. 2011), which allows 
us to rule out any simulated blends resulting in a com- 
bined color that is significantly redder or bluer than the 
target; (b) limits from the centroid motion analysis on 
the angular separation of companions that could produce 
the signal (Section 14. 3p : (c) brightness and angular sep- 
aration limits from high-resolution adaptive optics (Sec- 
tion 14. 1|) ; and (d) constraints from the measured tran- 
sit depth derived from our Spitzer observations, which 
place an upper limit on the mass (spectral type) of stars 
producing the blend. For the hierarchical triple sce- 
nario (in which the secondary star is physically bound 
to the Kepler targ et star), we considered dy namical sta- 
bility constraints (jHolman fc Wiegertl 119991 ). Surviving 
blend scenarios that satisfy all of the above criteria are 
folded into the "blend frequency": which is the prob- 
ability that such a finely-tuned blend lies near enough 
to the target star to be undetectable by adaptive optics 
imaging (using the sensitivity limits we find in Section 
14. 1[) . We compared this frequency with the expected fre- 
quency of true planets (the planet "prior" ) to derive the 
"odds ratio" . To estimate the planet prior, we employ 
the lis t of candidate planets (KOIs) from Batalh aTet al.l 
(|2013f) . restricted to main-sequence host stars. We as- 
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sume that this hst is complete (i.e., that all signals have 
been detected) and that the rate of false positives is neg- 
ligible (which assumption is justified by the findings of 
iMorton fc Joh nson 2011, who found a false-positive rate 
of < 5% for most KOls). 

We find that background eclipsing binaries comprising 
two stars furnish only poor fits to the Kepler photometry, 
and so are excluded from Kepler data alone. A portion of 
the hierarchical triple parameter space (with a star and 
larger planet gravitationally bound to the Kepler target 
star) provides good fits to the Kepler photometry, but 
these scenarios are then excluded by the combination of 
observational constraints described above. Background 
stars transited by larger planets, on the other hand, can 
mimic the Kepler photometry and remain consistent with 
the Kepler centroid motion, follow-up adaptive optics 
imaging, spectroscopy, and Spitzer constraint. We find 
that the frequency of background/foreground blends that 
satisfy these criteria is 4.19x10"^. The planet prior is 
estimated by counting the number of known KOIs that 
are in the same (Scr) radius range as the putative planet 
(105 in this case), and dividing by the total number 
of main-sequence Kepler targets observed during Q1-Q6 
(138,253). We obtained a planet prior of 105 / 138,253 
= 7.60 xlO~^. We conclude that a true transiting planet 
is 7.60x10"'' / 4.19x10"* = 18,000 times more likely 
than a blend, which allows us to validate Kepler-61 with 
a high degree of confidence. 

The simple procedure described above for estimating 
the planet prior does not take into account the period of 
the signal, which may be an important factor for small 
and long-period candidates (such as Kepler-61b) because 
such signals are rare (see below). Furthermore, the com- 
pleteness and purity of the KOI list of Batalha et al. 
(2012), on which our planet prior calculation relies, may 
decline with period and planetary radius, whereas we 
have assumed these concerns are negligible. These fac- 
tors may in principle influence both the planet prior and 
the blend frequencies we have just described (since we 
employ the KOI list not only to estimate the occurrence 
of authentic small planets, but also to estimate the rate 
of occurrence of larger planets in false-positive blend sce- 
narios). Therefore, instead of allowing eclipsing binaries 
and transiting planets with any orbital period to factor 
into the blend frequency calculation, we elected to redo 
the BLENDER analysis with a more realistic approach to 
allowed blends. First, we accept only blends with periods 
near the measured periods of Kepler-61 (within a factor 
of two) for both the blend and planet prior calculation. 
To address the concerns about completeness and purity 
of the KOI list, we performed separate Monte Carlo sim- 
ulations to establish incompleteness corrections for the 
KOI list and also to estimate the false positive rates for 
planets in the size ranges relevant to this calculation. A 
description of this work is forthcoming (Fressin et al., in 
preparation) . We obtained a revised frequency of back- 
ground/foregro und blends of 7 . 27 xl O~^. Examining the 
candidate list of lBatalha et al.l (|2013[ ). we found 22 KOIs 
in the relevant radius range with periods within a factor 
of two of the period of Kepler-61. Our simulations sug- 
gest that about 2.09 of these may be false positives, but 
also that the KOI list for signals of this size and pe- 
riod is in fact incomplete, requiring a correction factor of 



approximately 1.59 (i.e., a signal such as that of Kepler- 
61 could only have been detected around 63% of main- 
sequence Kepler targets). The corrected planet count is 
then (22 - 2.09) x 1.59 = 31.66. With this, the planet 
prior becomes 31.66 / 138,253 = 2.29 xlO"^. The final 
odds ratio for Kepler-61b is then 2.29x 10"'' / 7.27x 10"^ 
= 31,500, which 1.75 times more significant as we found 
with a more simplified approach. 

Blends that include a companion star <0.45 Mq would 
produce transit depths inconsistent with our Spitzer ob- 
servations, since they would produce transits depths 
more than 3cr deeper than we measure at 4.5 fim. These 
blends are thus excluded. For Kepler-61, the Spitzer 
results exclude all remaining possible physically bound 
configurations, which would have been the major cause 
of false positives otherwise. We depict this constraint in 
the bottom panel of Figure [11] 

We note that we have assumed that any signal with a 
SNR larger than 7.1 would have been recovered by the 
Kepler pipeline as a KOI, to compute the incomplete- 
ness correction factor. This optimistic hypothesis is a 
conservative one in our case, since a more realistic de- 
tection model would further increase the incompleteness 
correction to our planet prior. We conclude that Kepler- 
61 is an authentic 2.15 i?® planet with a high degree 
of confidence. We depict an illustration of the BLENDER 
constraints on false positives for Kepler-61 in Figure [TTJ 

5. DISCUSSION AND CONCLUSIONS 

5.1. Transit Times 

We depict the transit times of Kepler-61b in Figure 
[T2l We report no significant deviation from a linear 
ephemeris. 

5.2. Theoretical Composition of Kepler- 6 lb 
5.2.1. Bulk Composition and Atmosphere 

While we cannot estimate the mean density of Kepler- 
61b without a measurement of its mass, we can still dis- 
cuss plausible compositions, given its equilibrium tem- 
perature and radius. There now exist a sizable set of ex- 
oplanets with radii in the 1.0-3.0 R^ range with dynam- 
ically measured masses, though these span a large range 
of bulk densities from 0.7 g cm~^ in the case of Kepler-llf 
to 10.4 g cm~'^ in the case of CoRoT-7b (Kepler-68c may 
comprise a very dense exception, with p — 28_23i but the 
density range is large and relatively unconstraining). We 
list published masses, radii, and mean densities from the 
literature in Table 4. The two pl anets nearest to K epler- 
61 in radius are 55 Cancri e ( I Winn et al.l 120111) and 
Kepler-68b (Gilliland et al. 2013), the radi i of which lie 
within 0.15 R^ of the radius of Kepler-61. iCarter et al.l 
(|2012. ) found that, even within the small known sample of 
super-Earths with measured radii, a trend is apparent for 
planets with equilibrium temperatures < 1200 K: these 
tend to have "niini-Ncptune" compositions, with mean 
density <3.5 g cm^"^. However, these planet s (Ke pler- 
11b, d, e, f, & g, describ ed bviLissauer et al.| l201lL and 
GJ 1214b, described by iCharbonneau et alT i2009) are 
also near to or larger than 2 i?^, whereas all planets 
with measured densities higher than approximately 7 g 
cm"^ have radii smaller than 2 R^. It is therefore un- 
clear whether the low density of this set of cooler planets 
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is attributable to their planetary radius or their insola- 
ti qn, or both; t his qu estion is explored in greater detail 
in lWeiss et al.l (|2013[ ). who incorporated both quantities 
in relation to planetary mass in their exoplanetary mass- 
radius relation. 

We also consider th e theoretical atmosp heric content 
of a 2.15 Rq planet. iRogers et al.l ()2011|) consider two 
scenarios (core accretion and outgassing) by which plan- 
ets in the 2-4 Rq regime might retain a substantial 
hydrogen and helium envelope. Though that work fo- 
cuses on temperatures > 500 K, a cooler temperature 
would extend still longer the lifetime of a putative hydro- 
gen/helium envelope. For example, if Kepler-61 formed 
by core-nucleated accretion beyond the snow line, at 500 
K (substantially warmer), a hydrogen helium envelope 
fraction of 0.1% by mass is plausible for timescales <1 
Gyr. Alternatively, it the hydrogen content of the at- 
mosphere is outgassed from the planet (assumed to be 
formed from iron enstatite), a mass fraction of 1% by 
mass is plausible for timescales as long as 100 Gyr. These 
timescales (and their corresponding atmospheric mass 
fractions) should be considered lower bounds, given these 
formation scenarios, since a cooler planet like Kepler-61 
will retain an atmosphere for a longer duration, all else 
being equal. The synthetic planetary ra dius distribution 
.enera ted by the formation models of iMordasini et al.l 
2012f ). which assume a primordial hydrogen/helium en- 
velope, furnishes a good match to the Kepler candidates 
for planets larger than 2 i?^, but diverges from the Ke- 
pler results for smaller radii. This result may be at- 
tributable to the more terrestrial composition of planets 
smaller than 2 i?0 , for which the assum ption of a hydro- 
gen/helium envelope is no longer valid (jMordasini et al.l 
|2012). 

We conclude that a density larger than 7 g cm^'', which 
has only be observed for planets <2 R^, is unlikely for 
Kepler-61. Given its radius and comparatively low stellar 
insolation, its mass may be closer to the 3-6 g cm~^ den- 
sity range bracketed by 55 Cnc e or Kepler-68b wi th sim- 
ilar ra dii. If we apply the relation derived by Weis s et al.l 
(|2013D from the sample of exoplanets with radius mea- 
surements and masses <15OM0 (another power law ap- 
plies for more massive planets), we find a predicted mass 
and density for Kepler-61b of 3.2 M^ and 2.4 cm~^, 
respe ctively, near to that of Kepler- lib (iLissauer et al.l 

Mm- 

There also exist theoretical constraints on the sus- 
tainability of super-Earth atmospher es for higher mean 
molec ular weights. In particular, iHeng fc Kopparlal 
(|2012l) consider the stability of high mean molecular 
weight atmospheres belonging to super Earths orbiting 
low-mass stars in particular. The proximity of the habit- 
able zone to the star means that many super Earths will 
be spin-synchronized, with a permanent day and night 
side. In partic ular, the timescale for spi n synchroniza- 
tion is given bv iBodenheimer et al.l (12 0011 and stated in 
terms of orbital frequencv ft bv lHeng fc Kopparla (2012,): 



(7) 



_8Q_ /wN /Mp\ f a_ 
4517 \^) \M^) \Rp 



to lie within the range of 10-100 for rocky exoplanets, 
and in the 10^-10^ ran ge for gas giants, as stated in 
iGoldreich fc Soteilll966[ ) . Even with extremely rapid ini- 
tial rotational periods of the planet (e.g. 0.1 day) and 
values for Q which approach that of gas giants, Kepler-61 
is close enough to its host star where the spin synchro- 
nization timescale is less than 1 Myr. 

This poses a problem for atmospheric stability un- 
less the zonal winds' ability to redistribute heat in the 
atmosphere outstrips the radiative timescale of the at- 
mosphere. If this condition (namely, that the advec- 
tive timescale is shorter than the radiative timescale) 
does not hold, then the low temperature of "night 
side" of the planet can allow heavier elements to con- 
dense out, leaving the atmosphere unstable. An at- 
mosphere comprising heavier elements has a longer ad- 
vective timescale, since the wind speed is slowed as 
mean molecular weight increases (similarly to the sound 
speed). For this reason. Earth-like atmospheres (with 
mean molecular weights of 30) are particularly suscepti- 
ble to instability. The fact that Kepler-61 orbits a late 
K dwarf translates to a radiative longer than the advec- 
tive timescale (Hcng & Kopparla 2012), so an Earth-like 
atmosphere would remain stable. For later M stars (for 
example, an M 3.5V st ar, as adopted as a trial case by 
IHeng fc Kopparlall2012D . the radiative timescale at 0.25 
AU is shorter, so a 2.5 i?® planet possessing an atmo- 
sphere with mean molecular weight of 30 would be po- 
tentially unstable. 

5.3. Future Prospects 

We comment briefly on the feasibility of atmospheric 
characterization of Kepler-61b. It orbits a small star and 
may possess a hydrogen and helium atmosphere, both of 
which are favorable circumstances for transmission spec- 
troscopy. To perform a basic estimation of the expected 
change in transit depth at an optically thick wavelength, 
wc consider the atmosphere to be a ring with scale height 
-ff , where H = kTp/ fimg (and k is Boltzmann's constant, 
Tp is the temperature of the planet, /Zm is the mean 
molecular weight of the atmosphere, and g is the surface 
gravity of the planet). If we use a mass esti mate of 8 Mm 
for K epler-61b (near that measured by :Gil liland et al.l 
120131 for Kepler-68b, with a similar radius), then we ex- 
pect a surface gravity of 17 m s~^. If we assume the 
most optimistic case from a detectability standpoint, we 
also use molecular weight of 2 (corresponding to the 
hydrogen-rich scenario). Employing the the same equa- 
tion to estimate the change in transit depth attributable 
to the atmosphere as Miller- Ricci et al. (200^), we expect 
a change in transit depth given by: 



AD 



2TTRpH 
ttRI 



2RpH 



(8) 



where the planet's initial rotational frequency is given 
by w, Q is the tidal quality factor (and is believed 



which equates to AZ?=10 ppm, if we employ the values 
for Rp and i?* given in Table 2. If we instead assume a 
mass of 2.3 M^, , like that of the 2.6 Rq planet Kepler-llf 
(jLissauer et al.|[2011f) . then AD is correspondingly three 
times larger, at 30 ppm. This signal is approximately 
one-tenth the size of the 0.5 mmag values which might 
have been d e tecta ble in the atmosphere of GJ 1214b by 
iBerta et all ()2012D using the Wide Field Camera 3 on 
board the Hubble Space Telescope. However, Kepler-61 is 



12 



Table 4 

Properties of Transiting Planets from 1.4-3.0 ijgj with Dynamically-Measured Masses 



Name 


Radius 


Mass 


Mean Density 


Reference 






[R(b] 


[Me] 


[g cm-3] 






Kepler-68c 


0.953lHi2 


4 8 +2'5 
^■° -3.6 


9S+13 


Gilliland et al. 


(2013) 


Kepler-lOb 


1.486±0.035 


4 5fi+ll'' 
4 4c:+0.33 
4.45_o27 


S S+2-1 
"•o_2 9 

'•^"-0.59 


Batalha et al. 


r2011) 


Kepler-36b 


Carter et al. 


("2012) 


CoRoT-7b 


1.58±0.10 


7.42±1.21 


10.4±1.8=' 


Bruntt et al. 
Hatzes et al. 


(2010). 


Kepler-20b 


1 qi+0-i2 
-'-■^-'--0.21 


8.7±2.2 


"■^-2.6 


Gautier et al. 


f2012~) 


Kepler-llb 


1.97±0.19 
2.00±0.10 


4 3+2.2 
6.9±3.4 


3 1+2-1 
■^■1-1.5 

4.9±2.4 


Lissauer et al. 


(2011) 


Kepler-18b 


Cochran et al. 


(2011) 


55 Cnc e 


2.00±0.14 


8.63±0.35 


5.9111 


Winn et al. 


(2011) 


Kepler-68b 


2.31±0.07 


8.3±2.3 


3.32±0.92 


Gilliland et al. 


(2013) 


Kepler- 1 If 


2.61±0.25 


2 3+2-2 


'J- '-0.4 


Lissauer et al. 


(2011) 


GJ 1214b 


2.678±0.13 


6.55±0.98 


1.87±0.4 


Charbonneau et al. 


(2009) 




3S measurements of CoRoT-7b 








^ Differing ma; 


(lOueloz et al.ll200aiPont et al.ll2011ll 


furnish 



different mean densities; we have stated the most recently published values. 
Mean density ealculated from stated Icr limits in radius and mass. 



also 25 times dimmer in K band than GJ 1214, rendering 
the detection of an atmosphere around Kepler- 61 out of 
the reach of current instruments. 

Similarly, the radial velocity amplitude of Kepler-61b 
is increased by the small mass of the host star. In this 
case, assuming again a mass of 8 M^ for the planet and 
a mass of 0.64 Af© for the star, the planet induces a 1.8 
m s~^ motion of its star. However, though measuring 
a radial velocity signature of several meters-per-second 
has been achieved for dozens of exoplanets, these are all 
around very nearby stars. Kepler-61b, with its Kepler 
magnitude Kp of 15.0, is probably too dim for such study 
with current instruments. However, gathering additional 
Kepler observations of Kepler-61b will be helpful, par- 
ticularly given the fact that it will be observed in short 
cadence mode for Quarter 12 onward. 

5.4. Conclusions 

We present the validation and characterization of 
Kepler-61b, a 2.15±0.13 Ri^ exoplanet with equilibrium 
temperature of 273±13 K, orbiting a late K dwarf. We 
determine that the planetary hypothesis for the tran- 
sit signature of Kepler-61b is 30,000 times more likely 
than the false positive hypothesis, folding together ev- 
idence from high-resolution imagery, the stellar colors, 
the centroid position of the star from the Kepler images, 
the depth of the transit in the 4.5 /xm bandpass from 
Spitzer, and from the detailed comparison of the Kepler 
photometry to theoretical light curves of both planetary 
transits and stellar blends. Our measurement of the ra- 
dius and temperature of the star Kepler-61 is based upon 
a weighted mean of the directly measured radii and tem- 
peratures of a subset of nearby stars with the same spec- 
tral type, which quantities we apply as priors in our char- 
acterization of the planet. We present K-hand spectra 
and newly derived metallicities for this set of four similar 
stars, as well as for Kepler-61. The application of this 
empirical method, as compared to characterization from 
K-hand spectra and stellar evolutionary models, ulti- 
mately increased the size and temperature of the planet 
by 10%. We consider plausible compositions for Kepler- 
61b from the set of planets with similar radii and dy- 



namically measured masses, as well as from mass-radius 
relationships for exoplanets. We conclude that the planet 
is likely slightly too large to be terrestrial in composition, 
and likely possesses a significant atmosphere. 
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us. We thank Philip Muirhcad, Kathcrine Hamren, Ev- 
erett Schlawin, Barbara Rojas-Ayala, Kevin Covey, and 
James Lloyd for gathering, reducing, and sharing the 
TripleSpec iiT-band spectrum of Kepler-61. We thank 
the Spitzer team at the Infrared Processing and Analysis 
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the tenth Discovery mission. Funding for this mission is 
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Figure 1. ii'-band spectra for Kepler-61 (KOI 1361) and for 
nearby stars of similar spectral type. Overplotted in gra y are the 
metal-sensitive regions published by ' Mann et a l.' I|2013') ; we use 
the metric described in that work to calculate fFc/H] for each star. 
The H2O-K2 temperature index iRoias-Avala et al...2012. ) is also 
stated for each star. 
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Figure 2. Toil: Quarters 1-4 of Kepfer observations of Kepler-61. 
Dashed lines denote tlie intervals between quarters, and observa- 
tions in transit are indicated in gray, with vertical lines above and 
below transit times. Bottom: The results of the discrete correlation 
function applied to this portion of the light curve. The strongest 
peak corresponds to 36 days, which variability is present by eye 
in the 300 days of observations depicted in the top panel. Dotted 
lines indicate the strongest period and its harmonics. 
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Figure 3. FAST spectra of GJ 820B (top), GJ 380 (center) and 
Kepler-61 (bottom). We have denoted the wavelength regions em- 
ployed by spectral typing software, "The Hammer" , developed by 
ICovey et ail ||2007l') . 
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Figure 4. Stellar effective temperatures versus radius for Kepler- 
61 and for nearby standard stars, with stellar evolutionary models 
from Dartmouth (Cotter et al. 2008) overplotted, from ages rang- 
ing from 1-14 Gyr, at values of [Fc/H] equal to -0.25, 0.0, and 0.25. 
The directly-measured radii and temperatures of GJ 380, GJ 338A 
& B, and GJ 820 B are drawn from Boyaiian et al. (2012). We 
employ a radius and effective temperature for Kepler-61 from the 
weighted mean of these values, with an error bar that encompasses 
the highest and lowest means in the sample. The estimates for 
Kepler-61 from iC-band spectra (Muirhead et al. 2012a) an d from 
the broadband KIC photometry (^Dressing fc Charbonneaull201.j ) 
are also depicted. 
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Figure 6. Top panel: Kepler-61 transit light curve as a function of 
phase for Quarters 1-11, gathered in long-cadence observing mode. 
Bottom panel: Kepler Ught curve for Quarters 12, 13, and 14, 
gathered in short-cadence observing mode and binned in intervals 
of 7 minutes. Overplotted in gray is the best transit model light 
curve, with parameters given in Table 2. The effect of integrating 
over the 29.5 minute exposure time corresponding to long-cadence 
observations is apparent in the different shapes of transit during 
ingress and egress. 
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Figure 7. [Top panel): Sample of Kepler candidate exoplan- 
ets with radii < 3 H^, orbiting stars with T^g <4200 K, and 
with equilibrium temperatures < 400 K. The dotted lines de- 
note 270 K (above which run away greenhouse effect occurs, per 
IKaltenegger He Sasselovl 120111) , and 2 Hq. {Bottom panel) Sam- 
ple of Kepler candidate exoplanets with radii < 3 -R®, orbiting 
stars with T^g <4200 K, and incident flux levels within 3 times 
the value received at the surface of the Earth. The dotted lines 
here denote the h abitab le zone for stars from 2700-4500 K, per 
jKopDar apu et aI7 (2013). The number labels adjacent to points 
depict the KOI number. The v alues for planet candidates calcu- 
lated bv lMuirhead et all Il2012bl ) orbiting stars with Tgg <3800 K 
are depicted in blue. For candidates orbiting hotter stars, for which 
the X-band method tends to underpredict temperature, we depict 
values instead from Burke et al. (2013) in red (this stellar char- 
acterization relies instead upon the comparison of the broadband 
colors to stellar models of Pinsonneault et al. 2012). Where these 
latter values have been revised by Dressing & Charbonneau (.2013 ) 
using the Dartmouth stellar models (Dottcr et al. 2008), we have 
used those radius and temperature values and depicted the KOI in 
gold. The radius and temperature for Kepl er-61 (KOI 1361) from 
ii'-band spectroscopy reported by Kluirhead et al. (2012tji) is given 
by the blue error bar, while the revised Itr contour from this work 
is shown in green. 
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Figure 8. Left: For each iteration of the MCMC chain, we calcu- 
late the temperature of the planet at apastron (blue), periastron 
(red) and the time-averaged equilibrium temperature (green). We 
have assumed an albedo of 0.3. Right: Histogram of time-averaged 
equilibrium temperatures for Kepler-61b. The long tail of tem- 
peratures higher than 350 K is contributed entirely from orbital 
eccentricities > 0.8. 
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Figure 9. Left: K band adaptive optics image of Kepler-61. An 
additional companion is located 0.94" away from the target star, 
and is 2.9 magnitudes fainter. Right: The sensitivity limits to 
additional point sources in the neighborhood of Kepler-61 as a 
function of radial distance from the primary target. The filled 
circles represent the K magnitude limits and each point represents 
a step in FWHM away from the primary target centroid peak. The 
dashed line underneath represents the ii'-band limits converted to 
Kepler magnitude limits if a star were to have a nominal Kp-K 
color, as described in the text. 
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Figure 10. Le/f; Transit of Kepler-61b gathered with Warm 
Spitzer at 4.5 ^m, binned by a factor of 16. The best- fit transit 
model with depth derived from the Spitzer observations is shown 
with a solid red line, while the Kepler transit model (with Spitzer 
4.5 fim channel limb darkening) is shown in green. The Spitzer and 
Kepler transit depths are in excellent agreement. The transit depth 
we can rule out with Scr confidence at 4.5 /im is shown by a dashed 
red line. Right: The results of a residual permutation analysis on 
the Spitzer transit of Kepler-61b. We detect the transit with Scr 
confidence, and the depth inferred from the Kepler light curve, in- 
dicated by a dotted line, lies within one standard deviation of the 
depth inferred from Spitzer. 
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Figure 11. BLENDER goodncss-of-fit contours for Kepler-61b 
corresponding to the three different scenarios that contribute to 
the overall blend frequency: background eclipsing binaries (top), 
background or foreground stars transited by a planet (middle), 
and physical companions transited by a planet (bottom). Solid 
colored contours depict the difference in y^ between an authentic 
transit model and a blend fit with those parameters. Only blends 
within the solid white contour acceptably match the Kepler light 
curve (3o" difference in x~ between the blend and transit model 
fit; see lFressin et al.ll2011l . while red, orange, and yellow contours 
correspond to blend models disfavored by the Kepler photometry 
by 4, 5, and 6a-, respectively. The axes in each panel represent two 
of the dimensions of parameter space for blends. For the top two 
diagrams the vertical axis depicts the distance modulus difference 
between the two stars, while the horizontal axis corresponds to the 
mass (spectral type) of the putative secondary star. In the bottom 
panel (hierarchical triple scenario), the vertical axis corresponds 
the size of the planet transiting the companion star. The cyan 
cross-hatched areas indicate regions of parameter space ruled out 
because the resulting Sloan r'-2MASS K color of the blend is either 
too red (left) or too blue (right) compared to the measured color, 
by more than 3cr (0.15 mag). The green hatched regions indicate 
blends that are ruled out because the additional star is less than 
1 magnitude fainter than the target and would have been detected 
spectroscopically. Finally, the gray areas on the left represents 
the constraint from our Spiizer observations. The diagonal dashed 
green lines in the top two panels indicate the faintest blends that 
can mimic the transit: approximately AKp = 4.5 mag both for 
background eclipsing binaries and for background/foreground stars 
transited by a planet. 
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Figure 12. Kepier transit times for Kepler-61 from Quarters 1-14, 
as compared to the best linear ephemeris model. 



